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摘 X. 由 于 沙 尘 气 溶胶 散射 和 吸收 机 制 的 复杂 性 ,目前 还 没有 成 熟 的 可 应 用 于 沙 尘 天 气 的 大 气 长 
波 辐射 反 演算 法 。 通 过 对 比分 析 常 见 的 沙 侍 、 城 市 .乡村 、 海 洋气 溶胶 对 大 气 长 波 辐射 强迫 和 MO- 
DIS 通 道 辐 亮度 影响 的 基础 上 ,提出 利用 气 溶胶 光学 参数 对 线性 模型 进行 修正 ,构建 了 适用 于 沙 尘 
气 溶胶 条 件 下 的 大 气 长 波 辐射 估算 模型 ,并 利用 黑河 流域 4 个 观测 站 点 ( 花 寨子 荒漠 站 、 混 合 林 站 、 
黑河 遥感 站 和 张掖 温 地 站 ) 实 测 数据 对 模型 的 应 用 能 力 进行 检验 。 结 果 表 明 :4 个 站 点 遥感 反 演 的 
大 气 长 波 辐 射 均 方 误差 为 17.1~20.4 Wm? 18 2 H-12.3~-18 Wem, BT E T b b RYE GS 
光学 厚度 变化 和 辐射 强 迪 效应 ,修正 后 的 模型 可 显著 提高 沙 尘 气 溶胶 影响 下 大 气 下 行 辐射 的 反 演 


EX 键 词 : 大 气 长 波 辐射 ; 
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作为 地 - 气 系统 能 量 交 换 的 重要 表征 因子 , 近 
地 表 大 气 长 波 辐射 (Downward surface longwave radi- 
ation,DSLR) 在 能 量 收 支 和 气象 气候 变化 研究 中 不 
可 或 缺 。DSLR 是 整个 垂直 柱 中 大 气 散 射 .吸收 和 
发 射 的 结果 ,可 以 通过 复杂 的 辐射 传输 模型 (如 
MODTRAN) 利 用 大 气 成 分 垂直 分 布 状态 TT BE BG 
线 气 洲 胶 和 云 参 数 等 进行 估算 。 但 在 实际 使 用 
中 ,由 于 辐射 传输 模型 的 输入 参数 较 多 旦 难以 获 
取 , 容 易 引 起 误差 传递 。 且 由 于 DSLR 一 半 以 上 来 
源 于 地 表 100 m 内 的 大 气 发 射 辐射 , 受 近 地 表 大 气 
温度 梯度 影响 很 大 。 受 制 于 仪器 测量 精度 和 局 地 地 
理 信息 状况 , 当 研究 者 通过 地 面 观测 资料 特别 是 大 
气温 压 数据 参数 化 大 气 比 辐射 率 来 进行 DSLR 估算 
时 “7, 模型 稳定 性 和 普 适 性 较 差 。 近 年 ,部 分 学 者 
提出 利用 卫星 反 演 的 时 空 连续 性 和 多 层次 性 的 地 表 
定量 物理 参数 替代 建 模 时 需要 的 部 分 参数 估算 
DSLR™。 也 有 学 者 基于 辐射 传输 模拟 数据 建立 大 
气 层 顶 (Top of the atmosphere,TOA) 辐 亮度 或 亮 温 
数据 与 DSLR 数据 集 ,利用 透 过 率 和 气压 的 关系 代 
数 式 作为 通道 权重 函数 撒 述 各 层 大 气 辆 射 对 长 波 
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精度 ,减少 长 波 辐射 应 用 中 的 不 确定 性 ,对 于 旱 区 地 表 能 量 收 支 研究 提供 了 重要 参考 。 
沙 尘 气 溶胶 ; 辐射 强迫 ; MODIS; 黑河 流域 


下 行 的 贡献 ,建立 DSLR 和 传感器 辐 亮 度数 据 之 间 
的 统计 关系 模型 *'""。 后 者 具有 一 定 的 物理 基础 ， 
且 简 单 快 捷 ,算法 具有 较 高 的 稳定 性 ,适合 遥感 数 
据 产 品 的 业务 化 运行 ,但 未 考虑 云 和 气 溶胶 类 型 对 
模型 的 影响 。 

沙 侍 气 溶胶 作为 重要 的 辐射 强迫 因子 , 占 对 流 
层 气 溶胶 总 量 的 半数 之 多 , 沙 尘 气 溶胶 影响 下 的 
长 波 辐射 估算 是 一 项 具有 挑战 性 的 工作 , 沙 侍 光学 
性 质 和 辐射 传输 机 理 的 复杂 性 是 主要 误差 来 源 。 
近年 已 有 大 量 学 者 就 沙 尘 气 溶胶 开展 研究 , 茹 建 波 
等 尝试 从 气 溶胶 总 的 后 向 散射 系数 中 提取 沙 侍 
ES BOTA St RL, UNA) EB) FERREA 
别 利用 FY-4A、Landsat8 , VIIRS 等 开展 星 载 传 感 需 
反 演 气 溶胶 光学 厚度 的 研究 , 李 丁 等 " 利用 S5P/ 
TROPOM 数据 对 中 国 东 部 气 溶胶 单 次 散射 反照 率 
进行 反 演 。 相 对 成 熟 的 气 溶胶 光学 参数 反 演 方法 
为 研究 其 辐射 效应 奠定 基础 ,Xin 等 ”利用 便携 式 
太阳 光度 计 测 量 了 中 国 西北 4 个 沙漠 半 沙 漠 地 区 的 
沙 尘 气 溶胶 光学 厚度 和 不 对 称 参数 , 发现 其 光学 特 
性 和 辐射 特性 有 强烈 的 季节 变化 , 气 溶 胶 的 TOA 48 
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射 强迫 在 3.9~12 Wo 左右; 田 等 等 “利用 地 基 观 
测 资 料 计算 分 析 了 沙 侍 气 深 胶 对 大 气 逆 辐射 的 影 
i, HX ETE DF 0.3 AY, KAI EKARA 
浊 度 的 增加 而 增加 EE 33 Ts A ER DC A E T 
相反 的 趋势 ;Sicard 等 ”基于 辐射 传输 模型 对 沙 尘 


绕 该 区 域 已 开展 大 量 的 科学 研究 实验 ,积累 了 
大 量 水 文 气象 观测 网 数据 , 既 为 流域 生态 水 文 研究 
葛 定 了 很 好 的 科学 数据 基础 , 亦 可 为 本 文 研究 沙 尘 
气 溶胶 模型 提供 帮助 。 人 研究 区 范围 和 站 点 分 布 如 
图 1 所 示 。 


气 溶胶 的 长 波 辐射 强迫 进行 研究 ,发 现 沙 侍 粒子 的 
半径 .折射 率 、 分 布 状态 和 表面 温度 等 对 重 射 强迫 


1.2 数据 来 源 
12. 模拟 数据 ”考虑 到 大 气温 湿度 、 风 以 及 气 溶 


敏感 性 较 高 。 前 人 的 研究 集中 在 对 沙 侍 气 溶胶 参 
数 的 定量 提取 及 其 引起 的 辐射 效应 描述 和 分 析 上 ， 
对 沙 尘 气 溶胶 条 件 下 的 DSLR 遥感 估算 研究 较 少 。 
为 此 ,本 文通 过 分 析 不 同 地 表 参 量 和 大 气 廓 线 
场景 下 沙 尘 气 溶胶 变动 对 DSLR 和 MODIS 通道 辐 亮 
度 的 影响 ,利用 辐射 传输 模拟 数据 建立 TOA 辆 亮度 
和 DSLR 的 统计 关系 模型 ,发 展 新 的 DSLR 遥感 反 演 
方法 ,最 终 利用 黑河 流域 地 面 实测 数据 对 模型 进行 
验证 和 普 适 性 评价 ,以 提高 沙 侍 气 溶胶 条 件 下 
DSLR 遥感 反 演 精度 ,并 为 相关 研究 提供 参考 。 


1 数据 与 方法 


1.1 研究 区 概况 

本 文选 取 黑 河流 域 作 为 研究 区 ,地 理 位 置 介 于 
98°~102°E、38°?~42°N 之 间 ,面积 约 1.429X 10? km; 
黑河 流域 上 中 、 下 游 区 域 涵盖 草地 农田 湿地 等 多 
种 地 表 类 型 ,是 我 国 陆 地 河流 研究 的 重要 基地 。 围 


花 寨 子 荒漠 站 
黑河 遥感 站 


4 
Cex 
站 十 黑河 流域 边界 
海拔 /m 


胶 参 量 等 因素 处 于 动态 变化 中 ,地面 测量 数据 与 卫 
星 同步 观测 数据 难以 达到 同步 ,本 文 利用 辐射 传输 
模型 MODTRAN 生产 的 数据 集 分 析 和 评估 沙 持 气 溶 
胶 对 DSLR 遥感 的 影响 。 其 中 ,大气 廓 线 来 源 于 热 
力学 初始 分 析 资 料 (Thermodynamic Initial Cuess Re- 
trieval ,TICR) 数 据 库 ,每 条 廓 线 含 40 层 温度 .湿度 、 
压力 .高程 等 信息 。 考 虑 到 晴空 和 沙 尘 爆 发 的 场景 
因素 ,以 相对 湿度 和 连续 湿度 作为 租 选 条 件 选 出 
946 条 上 晴空 大 气 廓 线 ,大气 水 汽 含量 (Water vapor 
content, WVC) 范 围 为 0.05~6.15 g+ cm”, JE EE Ze ih 
度 范围 为 235~314 ,基本 嘻 括 了 研究 区 内 四 季 常 
见 的 廊 线 类 型 ;发 射 率 来 自 ASTER 地 表 发 射 率 库 ， 
针对 黑河 流域 地 表 履 盖 类 型 特点 ,选择 7 种 主要 地 
表 和 覆盖 类 型 进行 模拟 。 气 溶胶 对 辐射 的 影响 取决 
于 其 空间 分 布 .自身 的 物理 化 学 性 质 ( 包 括 粒 子 尺 
BE . 谱 分 布 .化 学 成 分 等 ) 以 及 光学 特性 等 (如 光学 
厚度 、 单 次 散射 反照 率 ) , 沙 尘 气 溶胶 相应 的 信息 取 


图 1 研究 区 范围 和 站 点 分 布 


Fig. 1 Distribution map of study area and sites 
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自 云 和 和气 溶胶 粒子 光学 特性 软件 (Optical properties 
of aerosols and clouds, OPAC) HAA”, Ayia ae 
普 适 性 ,基于 底层 大 气温 度 对 地 表 温 度 按 一 定 区 间 
变动 ,变动 区 间 范 围 为 -5~15 K, 步 长 为 5 kK。 

1.2.2 通 感 数据 ”本文 利 用 Aqua 卫星 上 的 MODIS 
数据 产品 来 验证 模型 的 可 行 性 ,采样 时 间 点 为 2018 
年 1 月 至 2019 年 12 月 。 主 要 使 用 6 种 类 型 的 MO- 
DIS 产 品 , 分 别 为 : 星 上 辐 亮度 产品 MYD021KM、 儿 
何 定 位 产品 MYD03、 水 汽 数据 产品 MYD05_L2. 云 
fits li = m MYD35_1L2、 地 表 温 度 和 发 射 率 产 品 
MYD11_12 以 及 气 溶胶 产品 MYD04(C6.1)。 

MYD021KM 是 1 级 辐 亮 度 产 品 ,包含 经 过 辐射 
校正 大 气 顶 部 的 反射 率 和 辐射 亮度 数据 。MYD03 
是 地 理 定位 的 1 级 产品 ,内 含 每 个 1km 分 辨 紊 像 元 
的 地 理 坐 标 等 信息 ,大 地 经 纬度 以 及 太阳 天 顶 角 和 
观测 角 、 卫 星 观 测 天 顶 角 (View zenith angle, VZA ) 和 
方位 角 海陆 掩 膜 等 重要 信息 。 

MYD05_L2 是 大 气 可 降水 量 2 级 产品 , 它 是 在 
表面 反射 率 线性 光谱 变化 的 假设 下 ,基于 多 个 不 同 
吸收 通道 反 演 水 汽 量 的 加 权 平 均 , 其 中 每 个 通道 的 
值 通过 查找 表 法 得 到 。MYD35_L2 是 大 气 2 级 标准 
数据 产品 ,内容 为 云 掩 膜 , 可 以 用 于 云 检测 并 进行 
睛 天 选取 。MYD11_L2 是 陆地 2 级 标准 数据 产品 ， 
其 中 包含 1 km 地 表 温 度 和 发 射 率 值 。 

MYD04(C6.1) 是 MODIS Collection 006 版 气 溶胶 
产品 ,该 产品 融合 深蓝 算法 和 暗 目 标 算法 的 合成 结 
果 ,并 在 气体 吸收 计算 、 瑞 利 散 射 计 算 .陆地 /海洋 判 
识 云 检测 和 产品 质 挖 方面 做 了 较 多 改进 ” ,空间 
分 辩 率 为 3 km。 

为 增强 模型 验证 的 可 信和 度 , 需 要 在 时 间 和 空间 
上 将 MODIS 观 测 数据 与 地 面 站 进行 匹配 。 空 间 上 ， 
为 保障 数据 的 有 效 性 ,匹配 过 程 中 站 点 对 应 像 元 周 
边 至 少 需要 有 10 个 有 效 像 元 。 时 间 上 ,地 面 站 点 
DSLR 数据 为 利用 卫星 过 境 前 后 15 min 的 测量 值 双 
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线性 内 插 得 到 卫星 过 境 时 刻 所 对 应 的 值 ,由 于 
DSLR 在 短 时 间 内 变化 不 大 ,该 值 近似 代表 观测 站 
点 所 在 MODIS 像 元 对 应 的 DSLR 值 。 为 了 消除 云 的 
影响 ,本文 只 选取 晴空 数据 。 为 确保 选择 的 像 元 达 
到 要 求 , 仅 选择 MODIS 温 度 产品 MYD11_12 质 量 控 
制 集中 标记 为 晴空 以 及 MYD35_L2 云 掩 膜 数据 判断 
为 晴空 的 像 元 ;另外 对 像 元 是 否 存在 气 溶胶 进行 判 
断 , 选 择 站 点 所 在 位 置 的 中 心 像 元 及 周边 存在 气 深 
胶 的 像 元 作为 符合 要 求 的 目标 像 元 。 
1.2.3 站 点 数据 本 研究 选取 黑河 流域 2018 年 1 月 
至 2019 年 12 月 的 站 点 数据 对 模型 进行 检验 。 数 据 
来 源 于 国家 青藏 高 原 科 学 数据 中 心 提供 的 墨 河流 
域 生态 -水 文 过 程 综合 遥感 联合 试验 站 点 中。 以 花 
寨子 荒漠 站 混合 林 站 、 黑 河 遥 感 站 .张掖 湿地 站 4 
个 不 同 下 垫 面 气象 站 ( 表 1) 数 据 进行 验 证 。 这 些 数 
据 是 由 地 面 辐 射 计 每 隔 10 min 观测 一 次 所 得 ,利用 
卫星 过 境 前 后 15 min 的 测量 值 双 线 性 内 插 得 到 卫 
星 过 境 时 刻 所 对 应 的 DSLR 测量 值 , 该 值 近似 代表 
观测 站 点 所 在 MODIS 像 元 对 应 的 DSLR 值 。 
1.3 研究 方法 
1.3.1 近 地 表 大 气 长 波 辐 射线 性 模型 算法 Wang 
等 "基于 透 过 率 和 气压 的 关系 代数 式 作为 通道 权 
重 函 数 描述 各 层 大 气 辐射 对 长 波 下 行 的 贡献 ,建立 
了 大 气 长 波 辐射 估算 的 混合 模型 法 。 该 算法 机 理 
明确 ,能 弥补 物理 模型 中 由 于 辐射 传输 过 程 过 于 复 
林 或 大 气 输入 数据 难以 获取 所 带 来 的 缺陷 , 且 直 接 
利用 TOA 辐 亮 度 估算 DSLR ,算法 计算 效率 高 ,稳定 
性 和 和 鲁 棒 性 较 好 ,适合 业务 化 运行 ”"。 该 模型 形式 
如 下 : 

DSLR, =a, *a,L,; + a, Lag + a,b) ta,L;, + 

a,L,, talog t aL 

式 中 :DSLRe 为 晴空 条 件 下 大 气 下 行 长 波 辐 射 估算 
值 ;ao~a 为 依赖 于 VZA 、WVC 以 及 地 表 高 程 的 回归 
ABs Lo Lo Hl Lails 9 MODIS 热 红 外 通道 测 得 的 


(1) 


#1 气象 站 基本 信息 


Tab.1 Basic information of meteorological stations 


气象 站 经 纬度 海拔 高 度 /m 下 热 面 类 型 [分 量 辐射 仪 位 置 
花 寨 子 荒漠 站 100.3201°E , 38.7659°N 1731 山 前 荒漠 6m 处 ,朝向 正 南 
混合 林 站 100.1335°E,41.9903°N 874 怪 柳 与 胡杨 24 m 处 ,朝向 正 南 
黑河 遥感 站 100.4765°E ,38.8270°N 1560 人 工 草地 1.5m 处 ,朝向 正 南 
张掖 湿地 站 100.4464°E ,38.9751°N 1460 芦苇 湿地 6m 处 ,朝向 正 南 
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TOA 辐 亮度 。 此 模型 可 以 解释 大 气 长 波 辐射 数据 
92% 以 上 的 变化 ,在 晴空 MODIS 的 数据 应 用 中 , 模 
型 的 标准 误差 约 16 W.m?, 但 其 主要 应 用 于 稳定 气 
溶胶 类 型 且 气 溶胶 含量 较 小 的 晴空 大 气 环境 中 ,在 
沙 尘 气 洲 胶 天 气 环境 中 的 应 用 能 力 还 未 得 到 检 
验 。 本 研究 拟 以 此 模型 为 基础 ,增加 相应 的 气 溶胶 
参数 ,拓展 模型 的 应 用 能 
1.3.2 算法 修正 思路 WERA F DSLR 的 反 演 难 
点 在 于 难以 量化 气 溶胶 对 大 气 的 长 波 辐射 强迫 和 
由 于 气 溶胶 变化 造成 的 通道 TOA 辐 亮 度 的 变动 。 
为 研究 沙 尘 气 溶 胶 的 长 波 辐射 效应 ,本 研究 首先 分 
析 了 4 种 不 同 气 溶胶 类 型 引起 的 辐射 强迫 效应 ,在 
此 基础 上 利用 气 深 胶 光 学 厚度 (Aerosol optical 
depth,AOD) 对 MODIS 第 29 .31.32 通 道 进 行星 上 辐 
亮度 校正 .优化 系数 ,构建 沙 尘 气 溶 胶 天 气 下 近 地 
表 大 气 长 波 辐 射 遥感 模型 。 

图 2 为 沙 尘 .乡村 城市 和 海洋 4 种 气 溶 胶 存 在 
时 的 气 溶胶 长 波 辐射 强迫 (Aerosol longwave radia- 
tive forcing, ALRF ) 均 值 偏 差 (Bias) 和 标准 差 (Stan- 
dard deviation , STD ) 随 AOD 变动 的 关系 图 。 从 图 中 
可 以 看 出 ,4 种 气 溶胶 的 ALRF 及 其 标准 差 随 AOD 
的 增 大 而 增 大 ,ALRF 与 AOD 之 间 近 似 呈 线性 关 
系 。 其 中 ,海洋 气 溶 胶 的 ALRF 最 小 ,在 AOD=1.0 时 
平均 为 3.85 W ^m^, KE STD X 1.1 Wem”; WEAR 


40r 
spot PRAN 
A 乡村 气 溶胶 
30+ * 城市 气 溶胶 
海洋 气 溶胶 
, 25 
Eo 
E 
10+ 
TE i 
中 到 SE S à A s i 


-5 
0.0 0.1 02 03 04 05 06 0.7 08 09 1.0 
AOD 


注 :A0D 为 气 溶胶 光学 厚度 ;ALRF 为 气 溶胶 长 波 辐射 强迫 ; 
Bias 为 均值 偏差 ,实心 符号 表示 ;STD 为 标准 差 ,柱状 符号 
一 半 表 示 。 下 同 。 
图 2 在 不 同 AOD 情况 下 4 种 气 溶胶 类 型 的 
ALRF 均 值 及 STD 
Fig. 2 Average ALRF and STD of four different aerosol 
types varies with the AOD 


胶 造 成 的 ALRF 最 大 , 随 AOD 变化 全 部 为 正 强 迫 ， 
ALRF 均 值 从 2.09 W -m° Zl] 29.8 W -m° , H STD 变化 
为 0.7~8.3 W.m。 

沙 尘 气 溶胶 条 件 下 的 MODIS 热 红 外 TOA 通道 
辐 亮 度 变 化 情况 如 图 3 所 示 ,TOA 通道 辐 亮 度 变 化 
ZHY Bias STD 和 AOD 同样 呈 线 性 关系 。 整 体 来 说 ， 
TOA 辐 亮 度 变化 率 Bias 较 小 ,但 通道 29、31 和 32 这 3 
个 地 表 温 度 通道 的 STD 相对 较 大 ,分 别 达到 3.51%、 
3.14% 和 3.32% ,说 明 在 改进 DSLRe 模 型 ,使 之 能 适 
用 于 沙 尘 气 溶胶 条 件 下 的 DSLR 估算 时 ,利用 AOD 


进行 辐 亮度 校正 并 进行 系数 率 定 的 思路 是 可 行 的 。 
6r 
4r 
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BT 01 02 03 04 05 06 07 T 09 L0 
AOD 
TE: TO A 为 大 气 层 顶 ; 实 心 符号 为 TOA 辆 亮度 变化 率 的 Bias ,柱状 
符号 长 度 的 一 半 为 该 变量 的 STD。 
图 3 沙 尘 气 溶胶 条 件 下 不 同 AOD 情 况 下 TOA 通 道贺 亮度 
变化 率 及 STD 
Fig. 3 TOA channel radiance change rate and STD of 


different AOD conditions under dust aerosol conditions 


在 此 基础 上 ,本 文 利用 AOD Xf MODIS 第 29、 
31,32 通道 进行 星 上 辐 亮 度 校正 ,并 优化 系数 , 沙 
侍 气 溶胶 天 气 下 近 地 表 大 气 长 流 辐射 遥感 模型 
(DSLR,) 形 式 为 : 
DSLR, = aq +aLyt+asLstas(bo XAOD+1)L,,+ 
a,(b, X AOD + D)L4 +a,(b,x AOD+1) (2) 
Ly, + a bos t aL 
式 中 :ao~ay 与 bo~by 为 由 VZA AOD 和 WVC 的 回归 系 
数 。 为 使 算法 更 具 代 表 性 和 减少 输入 参数 本 身 带 
有 的 误差 ,采用 在 WwWVC 和 AOD 相 邻 区 间 设 置 部 分 
重 闪 区 间 并 分 组 细 化 的 方法 来 优化 系数 。 其 中 ， 
VZA 分 为 6 种 情况 :0、33.56° .44.42° 51.32° 、56.25°、 
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60.00°; WVC(g*cm?)47 12 4H: [0.0, 1.0], [0.5, 1.5], 
[1.0, 2.0] [1.5, 2.5] 、[2.0, 3.0] | 12:5, 3.5] 、[3.0， 
4.0] 、[3.5, 4,5] 、[4.0, 5.0] 、[4.5, 5.5] 、[5.0, 6,0] 、 
[5.5, 6.5]; AOD 分 3 组 :[0.0, 0.4] , [0.3, 0.7] [0.6， 
1.01, 依 上 分 组 ,采用 最 小 二 乘法 共计 回归 216 组 系 
数 ,算法 构建 的 总 体 流 程 如 图 4。 


2 结果 与 分 析 


2.4 敏感 性 分 析 

为 评估 输入 变量 的 不 确定 性 所 导致 的 DSLR 反 
演 误差 ,需要 对 模型 进行 敏感 性 分 析 。 模 型 误差 主 
Z HH MODIS 仪器 噪声 和 校准 误差 、AOD 的 不 确定 性 
和 模型 回 有 误差 所 决定 。 模 型 总 体 误 差 e(ADSLR) 
可 表示 为 : 


e(ADSLR)= Ja(Losa) +6(AOD) +6(alg) (3) 


(Loss) TH (4) 
A aop ACOD) (5) 
式 中 : 5(Losm)、5(AOD) 25a 7g Hifi os BE LAI AOD 
不 确定 性 导致 的 DSLR 误差 ;6(alg) 为 模型 本 身 误差 ; 
a, 为 第 ;通道 对 应 情况 下 的 系数 ;2 为 第 ;通道 噪声 
误差 ;0 为 偏 导 ; ADSLR,,, 为 适用 于 沙 尘 气 溶胶 条 
件 下 的 模型 与 适用 于 晴空 条 件 下 的 模型 差 值 ;AOD 
为 沙 尘 气 溶胶 光学 厚度 值 ; A(AOD) 为 AOD 的 不 确 
定性 误差 。 根 据 Remer 等 3 对 全 球 尺度 MODIS 气 


8(A0D)- 


筛选 出 的 TIGR. 
大 气 廊 线 数据 


MODIS 通 道 
响应 函数 
观测 天 项 角 
(VZA) 
气 溶胶 光 
学 参数 


MODTRAN 4.0 


模拟 数据 集 i 
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溶胶 产品 的 真实 性 检验 结果 ,MODIS 气 溶胶 光学 厚 
度 产 品 的 误差 可 表示 为 : 
A(AOD)=+0.05+0.15AOD (6) 

根据 MODIS 的 等 效 噪声 误差 在 通道 27.28 .33 
和 34 为 0.25K, 在 通道 29.31 和 32 为 0.05 ;对 于 
WVC, FHF MODIS 产品 中 其 精度 为 0.53 gem"? ， 
在 VZA=60" 观 测 情 况 下 本 文 得 到 模型 的 总 体 误差 
(图 5)。 

从 图 5 可 以 看 出 ,VZA=60° 观 测 情 况 下 ,由 通道 
Ha z RE LFI AOD 不 确定 性 及 其 造成 的 (Lus) 和 
6(AOD) 随 着 AOD 的 增 大 而 增加 ,但 6(alg) 随 着 AOD 
的 增 大 而 减 小 。 由 于 算法 本 身 误差 在 整个 误差 中 
占 比 大 , e(ADSLR) 随 着 AOD 增 大 出 现 先 减 小 后 增 
大 的 趋势 ,最 小 值 出 现在 AOD=0.4 处 。 当 AOD= 
1.0, e(ADSLR) 达到 10.7 W.m*。 在 沙 尘 气 溶胶 条 
件 下 ,精确 的 AOD 值 对 DSLR 反 演 非常 重要 。 由 于 
不 同 输入 参数 的 误差 可 能 相互 抵消 ,实际 误差 可 能 
小 于 图 5 中 的 理论 估计 值 。 从 其 他 角度 进行 的 研究 
也 得 出 了 类 似 的 结论 , 且 e(ADSLR) 随 着 观测 天 顶 
角 的 增 大 而 增加 。 当 VZA=0, DSLR 均 方 误差 最 大 
为 9.6 W-m^, 
2.2 检验 结果 分 析 

最 终 科 选 出 在 花 寨 子 荒漠 站 ` 混 合 林 站 .黑河 
遥感 站 .张掖 湿地 站 的 匹配 数据 为 28、33、28 、29 
组 。 花 寨子 荒漠 站 AOD 的 变化 范围 为 0.06~0.83， 
WVC 变 化 范围 为 0.12~2.84 g+ cm? ; 混合 林 站 AOD 


修正 的 算法 


站 点 实测 数据 修正 的 算法 
MODIS 过 境 时 模型 的 敏感 性 
的 同步 数据 分 析 及 检验 


UE: TIGR 为 热力 学 初始 分 析 资 料 ;MODTRAN 为 大 气 辐射 传输 模型 。 
图 4 修正 算法 反 演 DSLR 的 总 体 流 程 图 
Fig.4 Flow chart of DSLR inversion by using the improved model 
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注 :VZA 为 卫星 观测 天 项 角 ; (Lus) 和 5(AOD) 分 别 为 辐 亮度 民 和 


AOD 不 确定 性 导致 的 DSLR 误差 ;6(alg) 为 模型 本 身 误差 ; 
e(ADSLR) 为 模型 总 体 误 差 。 
图 5 VZA=60° 时 不 同 AOD 人 情况 下 的 误差 


Fig. 5 Errors under different AOD cases when VZA=60° 


(a) 花 寨子 荒漠 站 
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AOD 


AOD 


的 变化 范围 为 0.09~0.86, WVC 变化 范围 为 0.29~ 
3.09 gcm ;黑河 遥感 站 AOD 的 变化 范围 为 0.09~ 
0.89, WVC 变化 范围 为 0.33~3.02 gcm ;张掖 湿地 
站 AOD 的 变化 范围 为 0.09~0.87, WVC 变化 范围 为 
0.44~3.43 g- em '。 图 6 为 卫星 过 境 时 两 种 模型 
DSLR 反 演 值 与 实测 值 散 点 图 。 

在 花 寨 子 芝 漠 站 , 沙 持 气 溶胶 天 气 多 发 生 在 春 
季 ,DSLR, 模 型 估算 的 均 方 误差 (RMSE) 为 20.3 W- 
m ,偏差 (Bias) 为 -5.8 W * m^, ix 22 Y [8] 2g -47.5- 
31.9 Wem; HEIR A bE AIR KAS REE 
春季 , DSLR, AIEA RMSE X 19.5 W mm, Bias 
为 -7.8 Wm ,误差 范围 为 -45.3~32.9 W «m^; FER 
河 遥 感 站 WE ORBE CERE TERE, DSLR BR 
型 估算 的 RMSE 为 17.1 W.m Bias 为 -1.8 Wem”, 
误差 范围 为 -36.5~29.3 W ^m^ ; CESK AEH , 沙 尘 
气 涂 胶 天 气 多 发 生 在 春季 , DSLR, 模 型 估算 的 


(b) 混合 林 站 
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注 :DSLR 为 近 地 表 大 气 长 波 锅 射 ;DSLR, 为 适用 于 沙 尘 气 深 胶 条 件 下 大 气 下 行 长 波 向 射 模型 估算 值 ;DSLR.。 为 畏 空 条 件 下 大 气 下 行 长 波 锅 射 


模型 估算 值 ;Bias 为 均值 偏差 ;RMSE 为 均 方 根 误差 。 下 同 。 
图 6 沙 尘 气 溶胶 条 件 下 两 种 模型 的 DSLR 反 演 值 与 实测 值 散 点 图 


Fig. 6 Scatter plots of DSLR inversion values and measured values of the two models under the condition of dust aerosol 
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RMSE Jj 20.4 W * m^, Bias Y-12.3 W* m^ , ix 25 YG [i] 
为 -46.9~18.6W.m?。 为 了 更 好 地 评价 两 种 模型 在 
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F DSLR 是 瞬时 值 , 且 研究 是 在 假定 DSLR 短 时 间 
不 变化 的 前 提 下 进行 的 ,在 采样 点 时 间 前 后 有 风 或 


沙 尘 气 溶胶 天 气 下 的 应 用 能 力 ,模型 估算 误差 关于 
AOD 和 WVC 的 散 点 图 如 图 7 和 图 8 所 示 。 
DSLRc 模 型 实际 数据 验证 结果 比 模拟 数据 验证 
结果 误差 大 ,主要 因为 模型 未 考虑 到 ALRF 以 及 没 
有 对 沙 尘 气 洲 胶 情况 下 受 影响 的 星 上 辐 亮 度 值 进 
行 必 要 的 系数 调整 。 沙 侍 气 洲 胶 不 只 产生 较 强 的 
ALRF ,其 对 星 上 辐 亮度 影响 也 较 大 ,特别 是 对 通道 
29、31 和 32, 最 大 变动 幅度 达到 3.51%。 从 图 6 中 可 
以 得 出 改进 的 算法 在 轻 度 和 重度 沙 侍 气 溶胶 条 件 


云 也 会 使 DSLR 产 生 较 大 侦 差 。 另 外 ,有 研究 表明 ， 
土地 者 盖 类 型 特别 是 作物 生长 阶段 的 发 射 率 变动 
对 DSLR 反 演 精 度 也 有 显著 影响 ” ,但 在 实际 模拟 
中 无 法 考虑 这 种 因素 也 是 误差 较 大 的 一 个 因素 。 
沙 尘 气 溶胶 颗粒 表面 不 均匀 性 和 角度 效应 也 不 容 
忽视 ,而 在 模拟 过 程 中 我 们 假定 气 涂 胶 为 球形 ,但 
实际 气 浴 胶 散 射 特 性 为 非 球面 ,也 会 影响 验证 结 
A. 


除了 以 上 因素 外 ,AOD .WVC 和 TOA 辐 亮 度 的 


下 表现 效果 都 比较 理想 ,由 于 气 溶 胶 对 星 上 辆 亮度 


不 确定 性 也 需要 着 重 考虑 。AOD 具有 空间 异 质 性 ， 


影响 并 对 其 进行 了 误差 修正 ,使 得 模型 在 一 些 高 
水 汽 的 廊 线 中 也 能 表现 良好 。 

但 也 明显 看 到 ,模型 误差 在 实测 数据 中 比 在 模 

拟 数据 中 大 的 多 。 原 因 是 多 方面 的 :在 验证 过 程 中 

卫星 数据 是 整个 区 域 成 像 时 刻 的 观测 结果 ,而 地 面 

站 点 仅仅 是 一 个 点 上 的 观测 数据 ,这 个 站 点 的 测量 

存在 系统 误差 以 及 随机 误差 ;本文 佑 算 的 沙 人 尘 天 气 
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受 地 形 湿度 .温度 等 的 影响 。 在 垂直 分 布 上 ,AOD 
和 WVC 各 层面 上 的 散射 和 吸收 均 不 同 , 其 垂直 廓 线 
上 的 变动 是 造成 许多 气候 模型 预测 不 准确 的 重要 
因素 。 利 用 激光 雷达 探测 大 气 参 数 的 时 空 分 布 和 
描述 气 溶胶 特性 是 提高 模型 精度 可 行 的 方法 。 但 
这 几 个 站 点 没有 测 得 对 应 时 刻 的 垂直 廓 线 数据 ,这 
对 模型 的 验证 及 结果 分 析 有 一 定 影 响 。 尽 管 OPAC 
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注 :DSLR, 为 实测 DSLR。 下 同 。 


图 7 沙 尘 气 溶 胶 条 件 下 两 种 模型 的 DSLR 反 演 误差 与 AOD 散 点 图 


Fig. 7 Scatter plots of DSLR inversion error and AOD for two models under dust aerosol conditions 
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注 :WVC 为 大 气 水 汽 含量 。 
图 8 沙 尘 气 溶胶 条 件 下 两 种 模型 的 DSLR 反 演 误差 与 WVC 散 点 图 


Fig. 8 Scatter plots of DSLR inversion error and WVC for two models under dust aerosol conditions 


提供 了 沙 尘 气 涂 胶 的 光学 特性 ,但 这 与 站 点 中 实际 
的 气 溶胶 光学 特性 还 是 有 差异 的 ,比如 源 地 、 气 溶 
胶 成 分 、. 粒 径 分 布 及 其 时 空 变 化 规律 不 同 ” ,也 会 
对 模型 检验 的 结果 产生 影响 。 另 外 ,云天 像素 被 误 
认为 晴空 像素 会 使 TOA 辐射 亮度 大 幅度 降低 也 是 
本 文 不 可 忽略 的 问题 。 站 点 符合 要 求 的 数据 不 多 ， 
模型 实地 验证 的 精度 可 能 受 个 别 样本 点 的 干扰 。 


3 结论 


(1) 该 算法 建立 在 以 AOD 校 正 部 分 敏感 通道 并 
进行 系数 率 定 的 基础 上 ,综合 考虑 了 AOD WVC 和 
VZA 等 因素 ,在 模拟 数据 垂直 观测 中 可 将 e( ADSLR) 
控制 在 9.6 W'm?。 敏 感性 分 析 结 果 表 明 , e(ADSLR) 
随 着 AOD 增 大 出 现 先 减 小 后 增 大 的 趋势 ,最 小 值 出 
现在 AOD=0.4 处 。 

(2) 在 黑河 流域 4 个 站 点 的 验证 中 ,改进 的 模 
型 RMSE 在 17.1~20.4 W * m^, Bias 介 于 -12.3~-1.8 
W.m?, 基 本 满足 实际 应 用 需求 ,但 站 点 检验 整体 误 


差 比 模拟 数据 中 偏 大 。 

总 体 来 说 , 相 比 DSLR. 模 型 ,在 WVC 和 AOD 及 
观测 条 件 已 知 的 情况 下 ,由 于 考虑 到 了 ALRF ,本 文 
提出 的 模型 能 显著 提高 沙 尘 气 溶胶 条 件 下 的 DSLR 
估算 精度 。 然 而 ,限于 目前 的 研究 条 件 和 数据 ,一 
些 方面 有 竺 改进 ,尤其 是 气 涂 胶 垂 直 分 布 剖 面 对 
DSLR 反 演 的 影响 有 竺 进一步 人 研究。 未 来 的 研究 将 
考虑 气 溶 胶 廓 线 变 化 对 DSLR 的 影响 ,并 利用 激光 
雷达 和 太阳 光度 计数 据 对 气 溶胶 进行 分 析 研 究 , 研 
究 气 溶胶 的 垂直 分 布 廊 线 ,定量 分 析 廊 线 变 化 引起 
的 气 溶胶 消光 后 向 散射 比 变 化 及 其 带 来 的 消光 系 
数 误 差 对 DSLR 反 演 的 影响 ;增加 站 点 数据 并 开展 
有 关 气 溶胶 影响 下 DSLR 反 演 模型 的 普 适 性 研究 ; 
深入 开展 土地 覆盖 类 型 变化 对 DSLR 反 演 精度 影响 
的 研究 。 
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Estimation of downward surface longwave radiation in Heihe River Basin with 
remotely sensed data 


WANG Zichao', WANG Chunlei, MA Junjun’ 
(1. School of Artificial Intelligence, North China University of Science and Technology, Tangshan 063210, Hebei, China; 
2. Consulting & Research Center of Ministry of Natural Resources, Beijing 100100, China) 


Abstract: Downward surface longwave radiation (DSLR) is indispensable for research about surface radiation 
balance. DSLR is one of the components of the surface radiation budget. Most of the DSLR inversion methods 
with remotely sensed data are on the conditions of clear sky. This method has a limited application range and 
their performances can be severely degraded in complex weather conditions as the effects of different aerosol 
types and levels are rarely considered. However, the arid and semi-arid area in northwest China is a frequent area 
of sandstorms, especially in spring each year, millions of tons of dust can enter the atmosphere, which not only 
causes significant damage to human life and production activities but also affects climate and the environment 
due to the physical and chemical effects of dust aerosol. It is difficult to study the optical properties of dust 
aerosols at home and abroad. Due to the complex mechanism of aerosol scattering and absorption, there is no 
mature DSLR algorithm that can be applied to dust sky, and this research is a tentative exploration in this field. 
For this issue, DSLR and top of the atmosphere (TOA) channel radiances under seven main land cover types are 
simulated to establish a simulation database. Secondly, the model coefficients are grouped by view zenith angle 
(VZA), water vapor content (WVC), and aerosol optical depth (AOD) to form a coefficient lookup table. Finally, 
based on the analysis of the effects of various four aerosol types and levels on atmospheric longwave radiative 
forcing and MODIS channel radiances, an improvement of the linear model using aerosol optical parameters is 
proposed to establish a DSLR estimation model constructed to apply to dust aerosol conditions, and four sites in 
Heihe River Basin, Gansu Province, China were used to verify the application ability of the improved model. 
Overall, the root mean square error of DSLR between the retrieved value and the site value is from 17.1 W +m” to 
20.4 W-m”, and the bias is from -12.3 W - m” to -1.8 W-m’” for the four stations. Due to considering the dust 
aerosol optical depth and radiative forcing of dust aerosol, the improved model is able to significantly increase 
the estimating accuracy of DSLR under dust aerosol conditions and can reduce the uncertainty of DSLR in the 
applications of radiation. From the results, the inversion accuracy of this model can meet the requirements of 
instantaneous DSLR research and serve as a reference for subsequent research. However, more case studies and 
comprehensive analysis are needed before the model is used to estimate DSLR in dust sky actually. 


Key words: downward surface longwave radiation; dust aerosol; radiation forcing; MODIS; Heihe River Basin 


